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I. INTRODUCTION 


A. BACKGROUND 

In modern times it often occurs that during battle, friendly ships are threatened 
by enemy aircraft. This thesis focuses on a ship’s ability to defend itself against aircraft 
with a 40 mm gun system. The computer program presented in this paper is a 
modification of a simulation written in the A Programming Language(APL) by Capt. 
C.Keeling. [Ref. 1] Keeling’s program simulates engagements between air defense gun 
systems firing non-explosive projectiles at maneuvering aircraft. In his thesis, Keeling 
suggested further research on proximity fuzed and contact fuzed high explosive(HE) 
rounds. This paper expands Keeling’s simulation to include the HE round with both a 
contact fuze and a proximity fuze 

The simulation shows how changing gun parameters can affect a gun system’s 
probability of kill against an aircraft. These parameters include: target range, projectile 
size, ballistic dispersion, muzzle velocity, firing rate, radar error, linear fire control, type 
of fuze (a new feature) and aircraft probability of kill given a hit or warhead lethal 
radius (a new feature). The use of this program increases the information available on 
the effectiveness of gun systems with and without HE warheads. The ultimate 
measure of success in the design and performance of an air defense gun system 1s the 
ability of the system to shoot down low-flying hostile aircraft within the combat 


environment. 


By “OBJECTIVE 

The objective of this thesis is to continue the development of a computer model 
that will be helpful in making decisions regarding the design and operation of an air 
defense gun system. In this thesis, the program is used to gain insight into how each 
of the design parameters of the gun system affect the probability of kill and which of 
the parameters are most important in the gun system design for air defense. In 


particular, the influence of the HE round on gun effectiveness will be studied. 
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C. APPROACHING THE PROBLEM 

The gun system is located on a ship at the origin of an X, Y, and Z coordinate 
system. The playing area for the model is the quarter of a hemi-sphere entirely within 
the first quadrant of the X, Y and Z plane as shown in Figure 1.1. The size of the 


playing area may be adjusted by modifying the program’s source code. 
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Figure 1.1 Playing Area of The X, Y & Z Plane. 
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The model assumes one-on-one engagements, no suppressive fire from the aircraft, 
and no terrain effects. The target is not  necessarilv attacking the gun system and 


may Or may not maneuver while within the gun's effective range. 


D. ORGANIZATION OF THE MODEL 

The air defense system consists of an early warning system, a tracking radar, a 
weapon control system and a gun with one or more barrels. The model studied here is 
a twin 40 mm gun using high explosive rounds located onboard the ship at the origin 
of the playing area. As an aircraft enters into the radar detection envelope of the early 
warning radar, the aircraft is detected and identified as hostile by the ship. The 
tracking radar is given the targets position. This radar tracks the aircraft and 
continuously computes the target range, elevation, azimuth, and velocity. The gun fire 
control computer solves for the lead angle in azimuth and elevation based on the 
calculated target positional data assuming no target maneuvering. The gun is 
positioned at the lead angle and fires a burst of proximity or contact fuzed high 
explosive rounds. The projectiles fly out toward the moving target. The target 
continues on its randomly selected flight path. The miss distance and assumed ballistic 
dispersion are used in the model to determine the probability of aircraft kill for each 
round Proc. For the proximity fuzed warheads, the projectiles do not have to hit the 
target to cause damage. The probability of aircraft kill given a detonation, Py y,Is 
determined using the Carlton diffused Gaussian kill function . For the contact fuzed 
warheads, the projectiles must hit the target in order to cause damage. The probability 
the aircraft is killed given the hit on the target, Py iz; 1s represented by the two- 


dimensional Carlton hit function. 


1. Scenario 

The general approach in this model is the development of the kill probability 
for one gun system firing one burst at one randomly maneuvering target. The 
probabilitv of kill for each round is used to compute a burst kill probability. The gun 
svstem computes target flight path characteristics, predicts a time-dependent intercept 
point in space, positions the gun and fires the projectiles. For contact fuzed HE 
rounds, the target vulnerability is represented by a vulnerable area (A_) in a plane in 
space that is perpendicular to the slant range (R) between the gun and the center of 
the target. For the proximity fuzed warhead, the target vulnerability is represented by 


a lethal radius. 


2. Model Input 
Miplitseeare Tequircdmo: the gun, the {Wze, the target, and ranges of 
engagement. Gun inputs are muzzle velocity, projectile coefficient of drag, burst size, 
rate of fire, and angular ballistic dispersion. The fuze input is the probability of fuzing, 
which 1s independent of range. Aurcraft inputs are velocitv, altitude, acceleration 
performance capabilities, and vulnerability data based on the tvpe and size of 
projectile and fuze being used. A flight path generates, for several different types of 


aircraft, maneuvers using a Monte Carlo simulation. 
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Il. PROXIMITY FUZED WARHEAD 


A. BASIC THEORY 

A proximity fuze detects the target and causes detonation of the high explosive 
warhead in the vicinity of the target. A warhead usually contains a powerful but 
relatively insensitive high-explosive that can only be initiated by the heat and energy 
from the primary explosive in the fuze [Ref. 3]. Proximity fuzes accomplish their 
purpose through “influence sensing” with no contact between the projectile and the 
target. These fuzes are actuated by some characteristic feature of the target, such as a 
reflected radio signal, an induced magnetic field, an interrupted light wave, a pressure 
measurement. or an acoustical impulse. “Proximity sensing” results in detonation of 
the bursting charge in the vicinity of the target or targets. A direct hit is not necessary 
to disable an aircraft and achieve the desired effect. An aircraft may be damaged in 
anv part of the large volume occupied by fragments from the exploding round. Since 
several air targets mav exist simultaneously, the sensing equipment must be able to 
isolate a selected target. Proximity fuzing had its origins in England early in World 
War Ii. Proximity fuzing is represented as VT fuzing (a code name used during Word 


War [I to imply variable tume fuzing) (Ref. 2 p. 76]. 


B. TARGET SUSCEPTIBILITY 

The modeling and quantification of the individual events and elements in an air 
defense encounter is referred to as a susceptibility assessment. Susceptubuility is the 
inability of an aircraft to avoid being damaged in the pursuit of its mission, to its 
probability of being hit [Ref. 2] The level of susceptibility of an aircraft in encounter 


with a threat is dependent upon three major factors: 


e The scenario includes the physical environment in which the encounter takes 
place, the air defense deployment and activity, and the aircraft flight path and 
tactics, including any supporting forces. 


e The air defense characteristics, operations, and lethality. 


e The aircraft observables or detectable signatures, any countermeasures used, the 
aircraft performance capabilities. and anv self-protection armament are the 
important factors associated with the aircraft itself. 
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One of the most important susceptibility measures is the closest point of 
approach or miss distance of a projectile to a maneuvering target. 

In general, the smaller the miss distance, the more likely the aircraft will be hit. 
Whether or not the aircraft is hit depends upon the presented area of the aircraft. This 
study assumed the presented area Ay of the aircraft is based on the six cardinal sides. 
The most the projectile can see in the intercept plane is part of three sides. The threat 
aircraft used in the study has a presented area from the area of the top Atop (50 m2, 
equal to the bottom), the area of the side Nee (30 m 2 the same on each side), and 


the area of the front Ar, (5 m*, equal to the rear). 


C. MISS DISTANCE 
As shown in Figure 2.1, the CPA is the mean value for an assumed Gaussian 


distribution of the projectile miss distance and the ballistic dispersion is the round to 


round variance about that point. 


Boye ee : 





Figure 2.1 Ballistic Dispersion, Bias, and CPA. 


The miss distance form for the round 1s given by 


Be + My = CPA (eqn 2.1) 


where 


© (Hy, oH) is the mean or bias in the x and v direction. 


De) THE EUZE 

The function of the fuze is to initiate the detonation of a warhead at a time and 
place such that maximum damage will be inflicted upon the target. Since the warhead 
damage volume ts preset, the achievement of maximum damage depends on the time of 
fuze initiation [Ref. 3]. The fuze package consists of a safety-and-arming device to 
keep the weapon safe until it 1s deploved and clear of friendly forces, a detonator to 
initiate the HE charge detonation, a device that senses the presence of a target (known 
as the target detection device (TDD)), and a logic circuit that initiates detonation at 
the proper time. 

The functioning of a fuze can be represented by the probability P, which may or 
mav not be dependent upon the miss distance. For example, a proximity fuze with an 
80% probability of detonating the warhead within a miss distance range r, of 100 ft 


and a 0” probability beyond r, can be modeled in the form 


PAx.y) = 0.8 when (x7+ yo)!" < 100ft 
PAx.y) = 0 when (x7+ yejle > 100ft 


EF. TARGET VULNERABILITY 

Vulnerability refers to the inability of the aircraft to withstand one or more hits 
bv the damage mechanisms, to its vincibility, to its lability to serlous damage or 
destruction when hit by gun fire [Ref. 2]. The critical components on an aircraft are 
those components which, if either damaged or destroyed, would lead to an aircraft kill. 
A general procedure has been developed for determining the critical components, their 
possible damage or failure modes, and the effects of the component damage or failure 


upon the continued operation of the aircraft. The procedure consists of 


e A selection of the aircraft kill levels or categories to be considered. 
e <n assembly of the technical and functional description of the aircraft. 


e The determination of the aircraft components of the aircraft and their damage 
caused failure modes for the selected kill levels. 
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1, State Transition Matrix Method 

The study by Keeling used the state transition matrix method to describe the 
damage state of the aircraft after one or more hits. This method allows the 
consideration of redundancy in the critical components and accounts for the fact that a 
single shot cannot kill more than one redundant critical component. This method 
States that a sequence of independent events can be modeled as a Markov process. 
iiesstate Space of the Markov process includes kill of the aircraft, no kill of the 
aircraft, a state space for a kill of each of the critical redundant components, and 
combinations of kills of redundant critical components that eventually lead to a kill of 
Mie adircrait., The state of the target at the begining of the engagement is represented 
as the initial state vector S?, One transition represents the state of the target after one 
random hit on the aircraft, and S] am represents the state of the target after the jth ave 


The matrix that accounts for the transition is |TM|. Thus 


sit] = pu} ) (eqn 2.2) 


The state transition matrix method is not appropriate for the HE round. 


2. Vulnerability Model for Externally Detonating Warheads 
For the proximity fuzed HE rounds, the aircraft's vulnerability can be 


represented by a Py yp function of x and y proposed by Carlton 
yy ) 
Eee eXP(-X~ fy )texp(-¥7/tp°) (ecm 2.3) 
The scaling parameter Tg can be related to the warhead lethal radius I); such that when 


r=ry, Py p=0.5. Collocating the one-dimensional Carlton kill function through 0.5 


When r=r gives rg= 1.20 1) as shown in Figure 2.2. 


F. Pcs 
This simulation uses the circular Carlton kill function method as presented by 


Ball [Ref. 2 p. 320]. The equation for Pre 1s 
> iy ae 4 5 
PeSS= rene apa “)exp(-(Hpy” a)-(Hby” B)) (eqn 2.4) 


where 


' 
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LETHAL RADIUS 





Figure: 2,2 PE ID Function of X and Y. 


SCs Cop eem 7) 
a i 2 ZB 
e fis (264, + 1°) 
2 (64,:54,) is the ballistic dispersion of the round tn the x and y direction. 
° Pe 1s a scaling parameter of the warhead lethal radius r, such that r= 1.2r,. 


e  P- is the probability of fuzing (assumed constant). 

The calculation of the probability of kill for a variety of ranges 1s made simpler 
by using angular errors for the dispersion errors. From the Figure 2.3 it can be seen 
that @,, and @,» are the angles representing the angular standard deviation of the 
dispersion errors about the closest point of approach. Figure 2.3 shows that 


Olen odues (eqn 2.5) 


where 


oP Odx l= Sdx2 Seeds 
® G4, 1s the angular ballistic dispersion in the x direction. 


By the same theory in the y direction can be shown to be 


18 


Sdy == GyyR 
where Ady is the angular ballistic dispersion in the y direction. 





Figure 2.3 Effect of Angular Dispersion with Range. 


Furthermore assumption is made that the angular ballistic dispersion as shown in 


Figure 2.2 is the same in both direction, hence 
dx = Gay = @ 
The equation for the Carlton kill function model, Equation 2.3, simplifies to 
Press = (tg?Pel ¥) * exp(-(Mpy? + Bpy?)/(Y)) (eqn 2.6) 
KSS = (Pe ¥ PU(Hyy 7 My DY qn 2. 
where 
¢  yis(2a?R* + 1°) 


The final form for the probability of kill is obtained by using Equation 2.1, giving 


Press = (1)7P¢ y)*exp(-CPA2/ Y) (eqn 2a) 


For this study, r is taken as five meters. Thus, rg Is six meters. 
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G. PROBABILITY OF KILL GIVEN A BURST 

The probability of an aircraft kill for a burst of rounds is directly proportional to 
the probability there is a clear line of sight from the detecting element to the aircraft. 
Given that the aircraft has been detected, the probability that one or more propagators 
will be fired at the aircraft is detonated by Pr. Thus, the probability that the aircraft is 
killed by one shot in an encounter, PRUE, is given [Ref. 2. p. 320] by 


hie a De ee Press (eqn 2.8) 


¢ Pp 1s the probability that the aircraft has been detected from the start of a 
search up to the present time t. 


¢ Pp is the probability that a propagator will be fired at the aircraft. 
e This study assumes that Pp* Pe = I. 


The probability that the aircraft 1s killed by the sequence of N gun shots ts unity 
minus the product of the individual probabilities of survival for each shot. Thus, for 


the burst of N rounds Py-c¢ in Equation 2.8 is replaced with 


Penuem ol 1 _ We, (eqn 2.9) 
i=] 


where the subscript 1 denotes the 1th shot. 
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{l. CONTACT FUZED WARHEAD 


Contact fuzes detonate the HE warhead upon contact with the target. 
Detonation can be made to occur on the surface or inside the target, depending upon 


the design of the fuze. 


A. TARGET SUSCEPTIBILITY 

The presented area a of the aircraft is defined to be the area of the aircraft seen 
Prete projectile in a plane normal to the path of the projectile at intercept. The 
procedure for determining the presented area seen by the round is determined based 
upon the difference between the projectile’s velocity and the relative velocity of the 


projectile with respect to the aircraft. This procedure is described in detail in [Ref- 1]. 


B. TARGET VULNERABILITY 

For warheads that do not use a proximity fuze, such as small arms and HE 
warheads with contact fuzes, a hit on the aircraft must occur to cause damage. Thus, 
the kill function that defines the probability the target is killed due to a propagator 
becomes Py 7; for random hits, and integration is carried out over the extent of the 
aircraft. In this approach, the probability the aircraft is hit by the propagator, Pry, is 
computed and multiplied by the probability the aircraft is killed given the random hit 
on the aircraft, Pp yy. The computation of Py; is developed by Keeling and 
accomplished using the Carlton hit function. The aircraft vulnerabilitv is represented 
bv its vulnerable area en cembered attire aim point, and any hit on the vulnerable 
area causes a kill. The computation for the vulnerable area is the product of presented 
area Ap and the probability the aircraft is killed given the random hit on the aircraft 
ery (Ref. 2] 


7) Pk iH Ap (eqn Bel} 
This model assumes the P;- ;; value for the round is 0.25, and Ay is given as 


im 9.25” Ap (eqn 3.2) 


This simulation uses the Carton kill function for the Pegg as presented in (Ref. 





Fa. 
2 2 
A, —Tp —TH, 
P poo = exp] ———- +. (eqn 3.3) 
9) 2 2 2 2.1 Z Z 2 Z 
(226, +Xq) (2n0, +yq) (200, +X) (200, + ¥q) 


The product of x, and yg can be taken as (1)equal to the presented area, Ap, or (2) 
equal to the vulnerable area, Ay (see Ref.2, PP 316-315). This study assumes X9 = Yo 
and uses Ap for the product. The means and deviations in Equation 3.3 are 
determined in the same manners as was done for the proximity fuzed round in Chapter 


II. Thus, the final form for the probability of kill given a single shot becomes 


Prec= (.25)Ap (21a°R2 + Ap)*exp{-mCPA*;(2ma7R2+ Ap? 3.4 
Ksg= (25)Ap (274 p )eXP{- p)} (eqn 3.4) 


D. PROBABILITY OF KILL GIVEN A BURST 
When multiple shots are considered, the single hit Py iz, 1s used for each shot. 
To determine the probability of an aircraft kill for a burst of rounds, the same 


procedure uses for the proximity fuzed warhead is used. Thus, 


Peipure = 1 - II 1 rss, (eqn 3:5) 
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IV. COMPUTER PROGRAM 


Keeling’s program is modified by adding subroutines PROXFUZE, CONTFUZE, 
BURSTKILL, and by changing the appropriate data to a 40 mm gun. This computer 
program consists of a main program and utilizes fifteen subroutines. The program was 
Written in the A Programing Language (APL). The computer program listing is shown 


in APPENDIN A, and a flowchart of the program is depicted in Figure 4.1. 


A. MAIN PROGRAM 

The main program uses a DRIVER program for initiating simulation and probes 
the user for inputs. The program DRIVER calls the GUN subroutine for simulating 
gun fire solution. After the simulation has been initiated, subroutine GUN simulates 
the air defense engagement by utilizing the appropriate subroutines necessary to 
simulate the engagement of an aircraft by the on-board ship air defense gun. 
Subroutine GUN collects statistics on the probability of kill and where the actual slant 
range of intercept occurred. These data are then sent back to the main program 
DRIVER. 


B. SUBROUTINES 


1. Initiation 
Subroutine GUN simulates the air defense engagement by calling the 
appropriate subroutines necessary to simulate the engagement of the aircraft by the on- 
board ship air defense gun. 
2. Generation 
The generation of a situation for air defense, aircraft's position, the radar 
located position with radar errors, and flight path maneuver of aircraft 1s subroutines 
RANGECHEK, RADAR, and AIRCRAFT. 
3. Calculation 
Subroutines INTERCEPT, MISSDIST, and MISSVECTOR are used to 
determine the predicted intercept point for the median round in the burst and the 


closest point of approach for each round in the burst. 
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Figure 4.1 Flowchart of Computer Program.. 


4. Susceptibility 
Subroutines APDATA, PKVECTOR, CONTFUZE, and PROXFUZE are 


used in determining the susceptibility of the aircraft. APDATA determines the 


presented area of the aircraft that is used in the Carlton kill function model in 
subroutine CONTFUZE. PKVECTOR stores the data for each round in the burst. 
5. Vulnerability 
SUbGoUtnes EK VECTOR. CONTPLZE, PROXFUZE, and BURSTKILL are 
used in determining the vulnerabilitv of the aircraft. Subroutine PRVECTOR stores the 
Sitaeier eaeh round in the probability of kill, Subroutine CONTFUZE and 
PRONXFUZE are used in determining the probability of kill for a single shot whether or 
not the projectile hit the target. 
6. Random Number 
Subroutines UNIF, UNIRAND, and NORRAND are used to generate the 
random number required in the simulation. Subroutine UNIF generates uniform 
random: numbers in the interval(a, b). UNIRAND generates pseudo-uniform random 


numbers. NORRAND is the generation of normal random numbers. 


C. SUBROUTINE PROXFUZE 


Subroutine PROXFUZE is used in determining the probability of kill for a single 
shot of the HE proximity fuzed round. This vulnerability data is used by subroutine 
BURSTKILL to determine the probability of kill given burst. 


D. SUBROUTINE CONTFUZE 

Subroutine CONTFUZE ts used in determining the probability of kill for a single shot 
of the HE contact fuzed round. This vulnerability data is used in subroutine 
BURSTKILL to determine the probability of kill given burst. 


E. SUBROUTINE BURSTRILL 
Subroutine BURSTKILL determined the probability of kill for a given burst. 


This probability of kill data uses the vulnerability of the aircraft. 
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V. ANALYSIS AND CONCLUSIONS 


This program model was initiated with the realization that the probability of kill 
is One measure of effectiveness that has a broad base application in the field of air 
defense gun system evaluation. This program for determining gun system probability 
of kill has simulated several parameters. These parameters describe the gun, target 
environments, and engagement procedures whereby one gun system encounters one 
maneuvering target. The data used is given in APPENDIX D. This study uses a box 
plot to illustrate the effects of several important parameters. This graphical 
imformation gives a quick sense for the distribution of Py for the individual test 
values of each of the parameters. The vertical box contains the middle fifty percent of 
the distribution, and the individual dots show all the outliers. The mean value of the 
probability of kill at each of the test value settings are connected with a line. The result 
of each simulation ts listed with the Empirical Statistics table in APPENDIX E. 


A. TYPE OF FUZE 

The probability of kill of the aircraft in a one-on-one encounter with both the 
proximuty fuzed warhead and the contact fuzed warhead, Py ip. is shown in Figure 
5.1. The effects that a proximity and a contact fuzed warhead have on the lethality of 
the 40 mm gun using a HE rounds can be seen in this Figure. The proximity fuzed 
warhead shows a higher probability of kill given a burst than the contact fuzed 
warhead against a maneuvering target. Consequently, the remaining studies 


concentrate on the proximity fuzed warhead. 


B. BALLISTIC DISPERSION 

Figure 5.2 shows an example of how the 40 mm gun probability of kill changes 
with angular dispersion. The degree of dispersion can immediately be seen, as well as 
how many of the data points were considered as outliers. Since the dispersion settings 
for optimality is arbitrary in this range, the assumption is made that the dispersion 


angle for the gun ts 1.0 milliradians. 


C. RANGE EVALUATION 
The evaluation of the Py for different range values was conducted for the 


maneuvering target. The maneuvering target was engaged by one burst for each range 
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Figure 5.1 Graphical Comparison of Fuze on 40 MM gun Lethality. 


value. The gun produced the highest probability of kill at the 700 meters, and the 
effective range was 1300meters against the maneuvering target. The air defense for a 


Teal situation requires about 2000meters for successful engagement. 
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Figure 5.2 Effects of Ballistic Dispersion on 40 MM gun lethality. 


D. EFFECTS OF TARGET MANEUVERS 

Figure 5.4 shows the effects that a maneuvering and non -maneuvering target 
have on the lethality of the 40 mm gun. Against the non-maneuvering aircraft, the 
probability of kill given a burst is about twice that for the maneuvering target. That 
means that the lethality of the air defense gun requires a non linear fire control 


predictor based upon the tactical air combat maneuver( TACM). 


E. EFFECTS OF RADAR ERROR 

Figure 5.5 shows the effect of radar error on the 40 mm gun. In today’s combat 
environment, the radar accuracy is an important factor. That means a large 
improvement in lethality can be obtained by reducing the radar error. The 100 percent 
error represents the current state of the radar system, so greater or less than 100 


percent reflects great differences in lethality due to technological improvements to 


today s radar system. 
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EFFECTS OF RANGES ON 40 Mi GUN 





Figure 5.3 Effects of Range on 40 MM Gun Lethality. 


F. EFFECTS OF MUZZLE VELOCITY AND RATE OF FIRE 

Increased rate of fire and muzzle velocity can significantly improve the 
probabilitv of kill for the rounds in the burst, as shown in Figures 5.6 and 5.7. 
Increasing the rate of fire decreases the time between projectiles, which means there are 
more projectiles in the area of where the aircraft is at intercept. Thus, the probability 
of kill is increased for each round in the burst. Increasing muzzle velocity decreases 
the time required to get to the target, and thereby reduces the effects of target 


Mmaanetiver. 


G. CONCLUSIONS 

This study has been conducted from a theoretical and intuitive point of view with 
conclusions withheld in many cases pending validation. The execution of this program 
demonstrated what appear to be some clear improvements 1s gun lethality. The sample 
results in each of the graphs show the difference in the probability of kill values. They 
represent a convenient and practical base from which to study or perform parametric 
analyses on existing systems or on svstems still on the drawing board. The 40 mm gun 


with the proximity fuzed warhead considered here for the given scenario and target 
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Figure 5.4 Effects of Maneuver and Non-maneuver on 40 MM Gun lethality. 


vulnerability data has an effective range of about 1300 meters against a maneuvering 
target. This gun can be improved by improving the radar, the fire control computer, 


projectile muzzle velocity, gun svstem rate of fire. These improvements should be 


attempted in the order given. 
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Figure 5.5 Effects of Radar Error on 40 \{M Gun Lethality. 


Another air defense gun system study should be conducted to determine the 
synergistic effects of changing more than one parameter simultaneously. The model 
should be further developed to account for non linear fire control. This can be easily 


done bv changing the radar and intercept section of this program. 


31 


OB 


0.6 


0.4 


3 
: 
3 
° 
& 


02 





Figure 5.6 Effects of Muzzle Velocity on 40 MM Gun Lethality. 
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Figure 5.7 Effects of Firing Rate on 40 MM Gun Lethality. 
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0.8897 0.028000 0.8617 
[6] AAVMIN<10.1408 0.1408 0.0988 00.859200.042100 
Oeso92 6.0527 0 0G 817 
C7] AVMIN<« 10.2018 0.2018 0.131700.7938200.070100 
0.7982 0.0702 0:0 0002 ea 
C8] aAVMIN< 10.2628 0.2628 0.164600.737200.098200 
0.7372 0.0982 0010 CO s6a9 
C9] a 'AV DATA FOR MAX VELOCITY SAIS? AStAWVE CLG hice 
Elana hihi 10.2073 0.2073 0.137200.7/ 207 Go Ce 
79270. 0701 0000.7226 
Cra pAVMAX«1 Ono 0.23410.1500.7659 Oho 3gse ore 
.-/659 0.084100 0 0 6s 
E12) AVMAX+ 10.2018 0.2018 0.1728 00. /S8Z20no. 112 ae 
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Wy 362.0.1192 O70 6.6028 
C13] aAVMAX< 10.3165 0.3165 0.190200.683500.126200 
0.6835 0.1262000 0.5573 
[14] a 'WHAT IS THE MINIMUM VELOCITY?! 
(15] VMIN<50 
(16] >+(VBUSTS100)/STOP 
[17] a 'WHAT ae TRE MAXIMUM VELOCITY? ' 
[18] MAXV<VO+VMAX 
b194 ke (BUST VMIN )x C(AVMAX-AVMIN ) )+ (MAXV-VMIN ) )+AVMIN 
->> 
i Pree: X<0 
e221 GO: 


V MARKOV; DIMENSION; AYMATRIX:; IDMATRIX 3 TRANSMATRIX 
© YNEW;COUNT:PM3L3N:X 

sega HE OG OLE ALICE 3 Bo oA a ree een 

axxx*x«* BUILDS LETHALITY MATRIX AND TRANSITIONS x%xx*x 

axx«x IT TO GIVE PROBABILITY OF KILL RxKKKKKKKK KKK 

DKKKK KKK KKK KKK KK KKK KKK KK KK KKK KK KK KK KK KK KKK KK KKK KKK KKK 


pe VULNERABLE AREA DATA AS AN ARRAY! 
< 


L<oX1 
DIMENSION<Lx0.5 
AVMATRIXN« (DIMENSION DIMENSION )ox1 


Pal 
C2] 
ie | 
C4] 
iS, | 
C6] 
ez 
8] 
10] nyo 84 THE INITIAL STATE VECTOR 
a 

bee) Aa'WRAT TS JRE PROBA Bory OF HIT 
fied aPa<D 
C14] Ca 
L15] PM<1-PH 
[16] TDHATRIX- ( tDIMENSION )°.x(iDIMENSION ) ) 
par |) eS aN (DIMENSION, DIMENSION )p ((1DIMENSION )*2)) 
C18] TRANSMATRIX« ( (PMxIDMATRIX )+AVMATRIX ) 
ae aa eee NUMBER OF MATRIX TRANSITION! 
[21] GO: Y<«(DIMENSION , DIMENSION )pyY 
22 | UNEW<TRANSMATE XxXY 
C23] YNEW<+4 (QYNE 
C24] RD  BIMENSION )pYNEW 
eo) Y<QYNE 
[ ei COUNT <W- al 
ea) N<COUNT 
C28] OE Ae 

V CUMKILL 
DME LE CEeCCCCCCLCCLLCeCSCLSSSCoCSLOCoSSSCLSSSSCoS SSeS eS eee. ee 2 
C2] ax* GIVES CUMULATIVE PROBABILITY OF KILL FOR xxx 
C3] axx A TARGET THAT RAS BEEN ENGAGED WITH MORE xxxx 
CU] Axxx THAN ONE BURST xxxx«**x 
C5 OK K KKK KKK KKK KKK KKK RK KKK ROR OR IOI OOK RO KOKO KK OOK 
cS = » (pRANGES ) )pP) 
C3] 
bel <* 8a) 
C10] Q<1-q 

WeneAb UNIF Ns A3B 
MOREE CCeCCCCCCC£CCOCCCLCeL SOS SSS SS SSS SSeS CLS SoS oS. eee 2 ee 
C2] ax GENERATES N UNIFORM RANDOM NUMBERS xxxxxxxkxkx 
C3] axxx ON INTERVAL(A,B) KKKKKKKKKKK KKK KKK KK KKK 
Reece ceceeceecceceeeCeeeceeerceeeceeec eee ee. ££ e ee a 2 2 2 2 2 
[5] A<AB(1] 
C6] peatte: 
C7] R<A+(B-A)x (C#21474U83647 )x7?7Np2147483646 
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FWP 
Lit ws 


omen | ae | eee | ee | 


‘roe (men | em | fe {em em | ee || | | 
PPRPPOO~IMONEWNPR 
POP OLLI St 


L_t_jl_t 


V 
V R«N UNIRAND B 


Qo Kw ok Oo ok ok kK ok OK ok OK OK KK OK OK ok ok ok kK kk ok ok ok ok kK OK KOK ok ok ok ok kK kk ok KOK Ow 


Axx CENERATES | NV PSEUDO UNE ORM, AIDEN NUMBERS xxx 
ene ee eee 2 oo eT Rte ee 2 le eo 


eo =p2 1 Jae (N? 100000000. yx 2XxBL24 )#100000000 


V Z<N NORRAND P:S8;3:I3T 
Q & & ok Ook ok ok ok Ok OK kok Ok OK Ok ok ok kK Ok ok ok OK ok ok ok ok ok ok ok ok ok ok ok ok ok ok Ook Kk kk kK Kk 


axxx*x GENERATES N NORMAL RANDOM NUMBERS xxxxx 


Qk Kk ok ok kk kK kk kK kk OK kK ok ok ok kk Kk ok kk ok ok ok OK kk ok ok ok ok ok OK ok KK 


P<? 


Z<«(N)p0 
I<1 
F10:7<2 UNTRAND O25 05 
T<(2xT )-1 
S<(701)]*2)4+702)]*2 
>F10x1S21 
ZLIIJ<P eta t sen x(TC1)xC(C 2x@S )#S)*0.5) 
ae 1) 


CORRECTIONS MADE TO KEELING’S PROGRA™ IN; 


Subroutine ALRORAP TT (36, 37,055) 
Subroutine MISSV EGTORA) 
Subroutine MISS IS Gro. sie ) 
Subroutine PROBHIT (6) 


40 


APPENDIX B 
APL COMPUTER PROGRAM 


1. NEW PROGRAM 


V DRIVERS COUNT ¢i 
DIR IORI OR ROR ORO RRR ROR ROR RIOR HORROR ORO ROO IO Ke 
Ak KO DATA INPUT Fe ORR OOOO KOO 
DR ROR OR OIRO OR OIRO ROR ORO ORR ORR ORO ORR OR ORO ORO IORI ORK 
2 ie nae DATA POINTS DO YOU WANT PER RANGE?! 

aA a ae 

'WHAT IS THE GUN FIRING RATE IN ROUNDS PER MINUTE?! 

" PIRERATE 600 
Sans THE MUZZLE VELOCITY?! 
Pee ee PAE weno sol Z2E C ' 
ene aise ANGULAR DISPERSION IN RADIANS?! 

— 

A 'WHAT IS THE DRAG FACTOR B=CDxVO0x0.5?! 


Deo. 

ae 'WHAT ARE THE a YOU WANT TO LOOK AT?! 
RANCES+* 1900 1300 100 

A eee ARE THE AREAS” oF Pie WAC A Le | 


AFR 

ASTDEL30 

ATOP<50 

RL<5 

PK/#<0.25 

Qo kk kkk Ok ow ok Ok ok kk Ok Ok Ok kk Ok kk Ok OK OOK Ok OO OO Oe OO Ok OO OK OK OK OW 
ax*x*x PROGRAM EXECUTION AND DATA COLLECTION FAO OCI 
Aw P IS PROBABILITY OF KILL FOR KOO 

Ak R SLANT IS RANGE OF TARGET ALT INTERCEPT xxx 
|? 2. O22. 2.2.2. 2.2.2.2. 2.2.2, 2,8, 2,9, 2, 9,2, 2,9,9,9,9, 9,9, 8,0,9,0,9,0,8,9,8,9,8,9,.8,.8 9.8, 8,9,8, 8.0.4 
R<10 

P<10 

I<1 

FINO EL ORR 
LOOP: RANGE<RANGES(TI] 


GUN 
Ber, EKDATA 
R«R ,RANGEL 


> (TECOUNT )/LOOP 


° > 22.22. 2. OS. O22 22.2.2, 2.2.2. 2.8,.2,.9,2,.8,9,8,9,9,0, 2, 0,20,9,8, 0,9, 9,9,9,0,8, 2,8, 9,9,9,8,.8,8,8, 2,8, 9 3 


axx DATA CONTROL VECTOR USED FOR GRAPHS IN GRAFSTAT xx 


ik a a ee ee EERE RE RA 


PROO TOWN FON H 


CONAN WVOPOWDODNHANEWOPOWONAMNLEWNHrPOWONOMEWNHROLA «As An aA 


eed eemmetnammenaieents cence ae a | ae | ms | Wow S| EGS |S | CS | | | W)C | | | MD OY | MO) 


(| | | | | | | | |e | | eee ee | | eS | | | | | || | | | | | | oe | ee | Ge’ ee | gee’ oes | ee‘ | Gee | Ge) 


FEFEF FF FE FWWWWWWWWWWHOMWNNNWNNNNNPRPPPBRPBPLE 


Ae aes 

A7 : 2 10 i 

A \RANGES* 

A 

ace ((DATASI), CE RANGES))pRANGES 
eee QC) 


V GUN: I: VAC:BODY;:TOP:SIDE3;VRX:VRY;:VRZ 
[oD OOO OO. we. OOOO. 00.0. 8.2, @ 2. OO RPO. 0.20.0, 8,9, 0,0,.0,9,8,9,0,0,8,0,9, 9,9, 0,0,0,0! 
Axxkkxxwk* STMULATES AIR DEFENSE GUN BY CALLING xxxxx 
Axxxwxwxwke*x APPROPRIATE SUBROUTINES xxxxxxxx 


2.2. 22.2.2. 2.2.2.2. 2.22.2. 2.2.2, 2,9, 9, 8,05, 8. 8,2, 2,9,2, 0,9, 0,9,0,9,0,9,2,9,9,9,.0,9,0,9,9,8,9,0,8,8,5 


RANGEI<10 
I<0 


Lm | | | |e | ge | 


Onl Wore 
LL tt 


4] 


PKDATA< 
LOOP3: RANGECHEK 
RADAR 
AIRCRAFT 
INTERCEPT 

Die oe 
RANCEI<RANGET ,RTGT 
ee a eae »PK 
SEE ETI 


LL 
Le 
S 
ay 


V RANGECHEK;&303;D;MAXALT;MINALT 


D222. 2.2.2.2. 2, 0,8, 8,8,8,8,0,8,0, 9,0, 0,8,.0,0,0,8,8,0,0,0,0,08,8,0,8, 8,8, 08,0,0,0, 8.8.2. oO. e eee, ee 


Axxx* DETERMINES A/C POSITION AT TIME OF FIRING xx 


Qo kk ok Ow OK ok Kk Kk ok ok ok OK ok ok ok kk ok kk ok kk kk kk Ok kk kk kw KK Kw KK KK OK OK 


A 'WRAT eT nee MAXIMUM ALTITUDE IN METERS?! 


C1] 

28 

(37 

C4] 

5) MAXALTS4 

[61 a 'WHAT is THE MINIMUM ALTITUDE IN METERS?! 
[71 MINALT<9 

[8] R<RANGE 

[9] 0<01.5707963 UNIF 1 

[10] +(R>MAXALT)/LOOP1 

[11] ZP<(MINALT.R) UNIF 1 

[121 +LOOP2 

[13] LOOP1:2P<(MINALT ,MAXALT) UNIF 1 
[14] LOOP2:D<((Rx2 Nob RD) bx *0.5 
[15] XF<Dx(20 

Licue vy 


F<Dx(100 
V 


V RADAR;S2;THETAA;R;DELRANGE; DTHETA; THETAE 


CA HIRO RRO ORO ROR ORO IORI ae 
[2] axx*x*x DETERMINES RADAR LOCATED TARGET BASED xxx 

C3] Ax*xx*x*x*x ON ERRORS INDUCED BY RADAR xxxxxxxxx 

| A*& kk KK KKK kK KK KK OK OK KKK Kw KK OK OK OK OK OK KK OK OK OK OK OK OK OK KK OK OK OK OK KK OK OK OK OK OK 
[5] $2<0.0012+(2+#RANGE ) 

(6] DELRANGE<1NORRAND 06.2 

eg) a et ; 

C8] DTRETA<2 NORRAND(O, 

C9] ely le a 

C10] R«RANGE+DELRAN 

C11] THETAA< 19 CRC ( (XPa2 )4(YF«2))40. 549 

(12) THETAA<THET 

(13] THETAE< 10 (ZF +RANCE) 

aire) ae ae AE 

C15] YR<Rx DTHETA L2]+THETAE ) )x eee ae CIJ+ZHETAA 
[16] XR<«Rx 50 DTHETAL2)+7THETAE) )xX(10(DITHETALIIJ+7HETAA 
id BETES 10(DTHETAL2)+Z7HETAE 


V AIRCRAFT; VMIN; AXMAX ; AXMIN; AZMIN;AZMAX:D;0;VZ1 
-VZ23;MINALT: MAXALT: AXY 34213422301 


: aniiitas THE MAX VELOCITY OF A/C IN METERS/SEC? ! 
1% Bese © ACCELERATION IN Z DIRECTION IN M/SEC2?' 


es Ba BRI I II II ROR IO I II FOR IC ACC CRE kK ke IR eRe oe 
C2] axxx*x* A/C PERFORMANCE DATA INPUTS xxxxxxxkxx 

{3 DK kw KKK KKK KK KK Ow Kw kk OK OK KKK OK OK OK KK KK OK KK KK KK KK OK Kk Kw Kw Ok 
C4] aa, eles al THEMINALT THE A/C IS ALLOWED TO FLY?! 

is. MINALT<¢ 

[6] a 'WHAT TS Aas MAX ALT THE GUN WILL ENGAGE? '! 

eran eae 

(8] aA Waa? 78 THE MIN VELOCITY OF A/C IN METERS/SEC?' 

C9] Vth 

[10 

elt 

alee: 

C13 
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pees A/C ee ee IN Z DIRECTION IN M/SEC2?' 
a 'MAX ACCELERATION DN XOR Y OFTRECTEOGN PN M/SECZs 

AXY<0. UxXAZMAX 

Se ee Oe tees Geek we wen de wx 4 Re 

axxk*x DETERMINES A/C FLIGHT PATH MANEUVER xxxxx 

Ok kk KK AT TIME OF FIRING wxwkxwkwkkkk 

DO ORR OKO ORO Rk kk kk kok 

VA<( VMAX ,VMAX) UNIF 1 

0+ 06.2831853 UNIF 2 


Zits (MTHEER ZEN 
VZ1<( (MINALT-ZF)+(MAXALT-MINALT) )xVZ2 
VARI WALD Oey A 
Vax? )-(VZ%2))x0. 5 
VX< Be ye Oe aC 
VY<(Dx(100[1] 
AZ1<(MINALT-ZF )+( (MAXALT-MINALT )+(AZMINx0.6)) 
AZ2<0.6xAZMAX 
AZe(AZ2.AZ2) UNITE 
01<06.2831853 UNIF 1 
Hecht (e013 
WeSC abev@)ak 


IMM MMMM Mama ma oama aMmamanmanrmamarrmnsa mrs 
WWWWWWWWWNHNYNNDYMNNDNYDNYNP RRR R 
OnNONFEWHOrPOW OAM FWHYOPRPOWONMOMHL 
aes Wane es | WON | WON WO | | WON Cy WO | CW MO] | |W WW WM) WO 

Ss 

aN 

ot 

rah ERO 


AY 
VRZ<1 NORRAND 0.1,VZ 
VRX<1 NORRAND ae Or x 
ae NORRBAND(VY,(0.1,VY 


V INTERCEPT 3: F1ITN;BURSTIME;:XF1;:YF13:ZF1:RPROJ:N;FIN 
DP * kw kkk Kk kk kkk kkk kK kk kk KKK KK KKK Kk KK KK Ok Kk OK KKK KKK KKK 
ax* DETERMINES PREDICTED INTERCEPT POINT BASED xx 
ax*xx ON LINEAR FIRE CONTROL COMPUTATIONS xxxxxx 
ax*x*x AND THE NEWTON-RAPHSON TECHNIQUE xxxxxxx 
DR AeCPTRERITE+cox... 
DELT<1+(FIRERATE+60 ) 
ister > a (ETRERATEs 60) 
N<_ (BUR 2D) 
TO<Nx BEEP 
TAU<0 
LOOP1:XPI<XR+(VRXxTO )+(VRXXTAU ) 
YPL<YR+ (VEYxTO een 
ZPI<ZR+(VRZxTO )+(VRZxXTAU 
RTGT< ( Geen )+(YPI«*2 )+(ZP 
VOxTAU )*+( (BxTAU )+1 
FIN<RTGT-RPROJ 
F17N< ( (XP TX VEX )+ (YPIxVRY ) 
F1TN<F17TN-(V04(1+(BxT7ZAU)x 
> (RTGT<3000)/LOOP3 
AIRCRAFT 
LOOP3: > ( CLEEN )<0.. 0001)/LOOP2 
TAU<TAU-(FINtF 1TN ) 
>+LOOP1 
ae $20 


ei) Or 
+(ZPIXVRZ ) )+RICT 
2)) 


‘oY gums’ somes Seen | Se) Sn Sen {mn | en Sn Geen en en Yen Yes Y en Y en on | On | Qn | oe en Oe | 
NNNNNPRRRPRPPRPE PRD ONONOEWNP 
FEWNORODONOM FW OWI WItIttity 


LI dd 


V MISSVECTOR 


0 OOOO ORO OIRO OR ROI ORO IRR OIRO AOI AOI 
[2] axx CALCULATES DISTANCE FROM A/C TO MEAN THEORETICAL x 

[3] axx INTERCEPT POINT FOR EACH PROJECTILE 1k AO I Ik 

MRS Se CCCCCCCCO CeCe SSS ES SSeS SSeS ee SSS SSS SS See 2S SSS ee, 
C5] APN<1i0 

[6] VIMPN<10 

L7] CPAN<10 

[8J] N<1 | 

(9) TAUO< 0.1 

[10] LOOP:MISSDIST 

C11 J oe KK RK RO ORO ROO ORO ROR ROKK kkk ok 
[12] ax*x* COLLECTS STATISTICS ON BIS, PRESENTED AREA, xxx 

[13] axx AND BURST VELOCITY. xxxx 

[14] axx* USED TO DETERMINE PKSS AND PK/BURST — xxxxxxx 
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Qo OK KK OK OK OKO KK OKO OOK KKK oo KK KO 
CPAN<CPAN ,CPA 

APDATA 

APN<AP 

VINENeV HPN, VIMP 


N<N+ 
> (SBURST)/LOOP 
>0 


mw ww aie dete tt 
NNN PRR Re 
NP OWDONOM 
| a | Ga | | Cs | | ey | 


VS ab dabei sF2:F3:KA3: YA: ZA: FIIN: PORN 
DKK KKK KKK KK KK KK KOK KK KK KK KK KK KKK KK KKK OK KKK KKK OK KK KK 
ax*x*x* DETERMINES TIME AT WHICH CLOSEST POINT xx 
Ak kk OF APPROACH OCCURS xxxxxx 

DQ ** KKK Kk kK kK OK KKK KOKO OK KK KK KK KKK KK KK KK KK KKK KK KKK KK 
VPX« XPI*RTGT )xVO 
VPY<( YPI#RTCT )xVO 
VPZ<( ZPI+RTGT )xVO 
LOOP1:BT<+(BxTAU )+1 
BT2<BTx2 
XP< PeYxTAU) Br 


LP<¢ 
TO0< 





N-1)xDELT 
XA<XF+(VXxTO )+((AX4#2 )x(TOx*2) )+0((VX4+(AXxT0) )xTAU) 
+(AX (TAUx2 


* 
YA<YF+ VYxPO a AY#2)x(T0Ox*2 ))+((VY+(AYxTO ) )xTAU) 


+(AY (TA 
ZaZP+ (UZxT9 +f AZ*#2)x(T0*2) )+((VZ+(AZxTO ) )xTAU) 
+(AZ (TAUx?2 
XAAP<XA-XP 
YAP<YA-YP 
ZAP<ZA-ZP 
F11<XAP*x 
F12<YAP*x 
F13<ZAPx 
SE eee 


J 
J 
] 
J 
J 
] 
J 
4 VEY eae 
FATN< 

J 

i 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

] 

J 


/PYtBT2 )} 
Lee 


Rn 
© 





A<BITx3 
F21<((F11#XAP X+(2xXVPXXBtA) )xXXAP 
ml 2xVPYxBt+A))xYAP 
+ 


2xVPZt*B+tA))xZAP 
EOCE2 


F22< 
F23<¢ 
F2TN<F21+F22+F23 

7h (TAU -TAUO )<0.0001 
TAU0N<TAU 

TAUSTAU- (F1ATN+F2TN) 


al< 

8S CLLOOP2*1I)230)/LOOP2 
>LOOP1 

LOOP? : LE 

TAUO 
/evO4( ( (Bx Tay 
ee VOxTAU ) 


2 ) 


+1 
{tt axe TAU )+1) 


V APDATA 3s VXIMP;VYIMP 3 VZINP CAPT ae 


CD) IK I RO RO ROK OOOO ROKK IOI RRR IOI Ok 
[2] ax*x COMPUTES PRESENTED AREA OF THE AIRCRAFT xx 

ECR MEh See. ee eee ee ee eee oe 2 2 oe 2 2 eee eee Se eS eS See ee ee! 
Cu] VXIMP<VX- (VPxXP+RPROJ 

[5] VYIMP<VY-(VPxYP+RPROJ 

[6] VZIMP<VZ-(VPxZP+RPROJ 

ia) VINE ( (VAINP«3 )+(VYIMP2 )+ (VZIMPx2) 0. 5 

[8] >(VIMPSVP)/LOO 

[9] AP1<ASIDE- APTN - VP)x(ASIDE-AFRT )+# vary) 

[10] AP2<ATOP-( (VIMP-VP)x(ATOP-AFRT )+#(|(VA))) 

C11] +LOOP2 

EAD I He, SB: AP1<((VIMP-VP)x (ASIDE-AFRT)+#((|(VA))))+ASIDE 
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fSE ee EE Lee eee 


3] AP2<((VIMP-VP)x(ATOP-AFR 
i] L00 oe spetape apes OTE (EEG RG sO DD Dh 7-046) 2 


WIP RVECTOR: 7 3 NV 


MMe b.t.2.0.0,.0.0,0. 8.0 2 Ce eee eee ceeeLeeeCeeeCeCeeCeeeeeereee ee £2 2.23 
[2] axx STORES PROBABILITY OF KILL FOREACH ROUND xx 
Ct 82.0,0,0,0,.0,0,. 0.0.0.0. 0.0.02 2 2 eee eee eee ee oe fee ee te eee oe eee eee 
C4] PKSSN<10 
[5] N<pCPAN 
fen 2 <1 
m7 ZOOP:CPA<CPAN[T] 
See kt he tee ee CeCe Ce CSS SS SCCLLLSoLSSSSCCSoLoSoSoLoS eS SSeS ee eee 3 
C9] axxx*x OPTIONAL SELECTION FOR TYPE OF FUZE xxxkxxxkx«x 
L10] axx*x FOR PROXIMITY FUZE, USE PROXFUZE FORK Ak Ok a 
[11] axxx FOR CONTACT FUZE, USE CONTFUZE OOK 
eek x KK RK KK RRR IR RK RR KEK AR KAR KAKAK RAR KER KKK ARK RK 
mie) PROXFUZE 
[14] aCONTFUZE 
mie eee IA Gaye. 
[17] * (758) /LOOP 
V PROXFUZE;A3B3C3D3R3RL3RO3AL3 PF 
iA) O&K KKK Kk KKK KKK KK KK KK KK KK KK KK KK KK KKK KKK KKK KKK KKK KK KK 
Meeax* PROXIMITY FUZED WARHEAD STORES PROBABILITY xx 
C3] axxx OF KILL FOR SINGLE SHOT xxxxxkk kkk 
Cb] KK ROR ROR OR OR OIRO ROR ORO ROI I kk kok 
[5] R<RANGE 
Pel PF<1 
C7) RL<5 
C8] RO<(1.2xR 
[9] A<(CROx2)xP 
[10] (ROR2 OEE nun) )4(ROKD) 
me) C<At 
p12. Dis 1x (CPA*2))#((2x (DISP*2 )x (R*2))+(R0*2)) 


ee (xD) 


V CONTFUZE;SA383C;D3R 


] D<e(-3. Na NG ah ae W214 
J PESS<Cx (xD) 


Gn. 4.4 2.2,%,5,5,0,0,0,2,0.0,0,0,0, 0,0, 0,0, 228 228 eee 888228 e eee 8 eee eee 
[2] axxx CONTACT FUZED WARHEAD CALCULATES xxxxxx* 

C3] ax*x*x PROBABILITY OF KILL FOR SINGLE SHOT xxxxx 

C4] axx** AV=PK/HxXAP, ASSUMED PK/H=0.25 xxxxkkkx xk 

Mol Wet kt.4.2,%,4,4,5.0, 0.0, 0,0,0,5,0.5,5,5,0, 0,0, 0,0,0,0,0,8,2,0,2,0,8, 0,25, 5, 0,8, 8,8 2 2 eee! 
fel iE oe 

moe, A<0.25xAP 

C3] Be (3x3. 142 )x (DISPx2 )x (Rx2 ) )+(AP*2) 

ee C<A+B 

beg, 

[11 


V BURSTKILL:AK3:BK 
D&K KKK KK KK KKK KK KKK OK KK Kk KK KK kK KKK KKK KK KKK KK KKK KK KK 
ax*x*x PROBABILITY OF KILL GIVEN BURST xxxxkxkwexwewex 


Q ok ok kk ok ok kok kkk kok kok kkk kok kk kk OK KK kk KKK kk KKK Kk KK 


RAK< ( ee ODA LATIOVipeksS 
AK<PKSSN 


BK<1i-AKk 

Chex /EK 

EK<1-CK 

PK<EK 
ORO IORI OIRO ORI ROK AIK 
J] Axkkkx PE ae OF KILL ENCOUNTER xxx xxxxexx 
{ Qk*kKkKKK*x ASSUMED PDXPFHA1 RRR KKKKKKKKKKKKK KKK 


aAxk*x*k*« PK/E EQUAL TO PKBURST xkKKKKKKKKK KKK KK KK K 


Mmmm nmr 
Ree OONMnFWNrF 
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SYMBOL 


APPENDIX C 


SYMBOL NOTATION 


DEFINITION 


eNITS 


sie siecle che ote he oe che cic 2% Se tis oe oye ce 2s oe ce Se oe ob ae se ce ie 8% cl aye a ok 0% 9% 3s sis 9% oe oe cis 0% of oie 0% cle of of 9% fs oh ok ok 3b ok o& se oe a oe oie oe Se ee che ch oh oe Ob oe ok 5% ce 


area of aircraft front 
area of aircraft top 

area of aircraft side 

max aircraft acceleration 
min aircraft acceleration 


aircraft presented area 
aircraft vulnerable area 
PAmCeIOl ay 

acceleration of aircraft 
coefhicient of drag 
gravitational constant 

gravitv drop 

lethality matrix 

probability of aircraft hit 
probability of component hit 
probability of kill 
probability of kill engagement 
probability of kill given hit 
slant range 

projectile range 

time 

aircraft total velocity 

burst Velocity 


muzzle velocity 
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meters” 


2 
meters 

2 
meters 


2 
Meters, sec 
2 


meters, sec 
4 
meters~ 
y 
meters 
Rene. 
meters sec 
ees 
meters;sec 
1, meters 


= 
MIELE TS SEC = 


meters 


MELELS 
meters 

Sec 

Meters Sec 
meters, Sec 


IMETETS, SEC 


Vp 


Vx’ py’ pz 
Vo ene 


\ Vee 


84-95. 
Q..0 


projectile velocity 

projectile velocity component 
radar aircraft velocity 
aircraft velocity components 
Fanse Ofeve 

uncorrected wind 

aircraft position at time t 
aircraft position at firing 
projectile location 

radar located target 

angular ballistic dispersion 
projectile displacement 
range error 

azimuth and elevation error 
linear mean bias 

gun azimuth and elevation 
projectile linear range 
linear ballistic dispersion 
projectile flight time 
angular mean bias 


target azimuth and elevation 


4§ 


meters/sec 
meterssee 
meters! sec 
meters. sec 
meters sec 
metersesec 
meters 
meters 
meters 
meters 
radians 
meters 
meters 
radian 
meters 
milliradians 
meters 
meters 
sec 
radians 


mulliradians 


APPENDIX D 
EQUATIONS AND DATA 


I. GOVERNING EQUATIONS 


a. The Siacci Formula 


This equation used here to determine the velocity of the projectile at a point in time is 
V_ = Vo! (C,Vol 224 1)" D.1 
p ee) ( dO t ) (eqn ; ) 


Where these two equation can be seen [Ref. 1]. 


e Cg 1s the coefficient of drag of the projectile. 


. Vy is the velocity of the projectile. 


e =Vo is the muzzle velocity. 


e Tis the time of flight of the projectile. 


b. Range of the Projectile 


The formula for determine the range of the projectile at a point in time is taken as 


R, = Vg t(CgVol 22+ 1) (eqn D.2) 


The derivation of these formula can be found in [Ref. 4] 


c. Gravity Drop 


GD = (-1,6)g "(1 e2QVai Vg) 4) (eqn D.3) 
2. DATA 
This study looked at 40 mm gun system, and firing high explosive round with 


both a contact fuze and a proximity fuze. The unclassified svstem characteristics were 
obtained from [Ref. 5] 
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All references to vulnerable area data for specific air threats to these guns are classified, 
therefore aircraft probability of kill given a hit and warhead lethal radius is entirely 
svntheuc. Vulnerable area data for many different aircrafts and many different gun 
svstems is readily available for personnel with appropriate clearances. The optimal 
angular ballistic dispersion was determined using the computer simulation and the 
coefficient of drag for the projectiles are assumed,the other data was from Jane’s 
Weapon System. The aircraft's vulnerability measure r is assumed as discussed in 
Chapt 2, and Ay is determined by using the probability the aircraft is killed given the 
random hit on the aircraft Py 1; as discussed in Chapt 3. 


WADE 
40 MM GUN CHARACTERISTICS 


firing rate Ges / Mach | 600 
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APPENDIX E 
EMPIRICAL COMPARISON OF MARGINAL DISTRIBUTIONS 


Wee FYPE OF FUZE 
a. Proximity Fuzed Warhead 


TABLE 2 
PROXIMITY FLUZED WARHEAD 


PARAMETER VALUES: (EOC je 7 OCRE aeons (1900) 


POPULATIONS Roles Rolo RG 
WEIGHTS 0 
X LABEL : RANGE IN METERS 
Y LABEL : PROBABILITY OF KILL 
POPULATION NO. OF  # waweecceee PERCENTILES----------- 
NUMBER POINTS Y MEAN Y STD DEV. 0.25 0.5 es 
100 50 0.80613 0.1907 0.72854 8.6001E-1 9.4777E-1 
700 50 0.2179 0.23669 0.021002 1.3836E-1 3.2233E-1 
1300 50 0.051488 0.15441 0 9.2321E-13 7.8509E-3 
1900 50 0.0094727 0.065348 0O 0.0000E0 6.1062E-160 
POOLED STANDARD DEVIATION ESTIMATE: 0.17357 
b. Contact Fuzed Warhead 
TABLE 3 
CONTACT FUZED WARHEAD 
PARAMETER VALUES: (100);(700);(1300);:(1900) 
POPULATIONS : CR14:CR13;:CR12;CR11 
WEIGHTS : 0 
X LABEL ; RANGE IN METERS 
Y LABEL : PROBABILITY OF KILL 
POPULATION NO. OF | waeeeecce-- PERCENTILES---------- 
NUMBER POINTS Y MEAN Y STD DEV 0.25 0.5 0.75 
100 60 0.017203 0.0056854 0.013458 0.016796 1.9509E-2 


0.0092551 0.0053786 0.0055805 0.010485 Ne E72 
O7005507 3°70 CR O00GSS096 "5 32570555 
0 0 PP O420en-15 


700 60 
Je 16f0 60 O,00S tT 
1900 60 G©2C004S311 0.002209 


POOLED STANDARD DEVIATION ESTIMATE: 0.0049179 


Ze BALLISTIC DISPERSION 


TABLE 4 


BALLISTIC DISPERSIG 


PARAMETER VALUES : (0); (0200057770. ee (0.0015) (0. doz 


POPULATIONS AP1;AP2;AP3; -AP4;:A 
WEIGHTS : O 
X LABEL : ANGULAR DISPERSION IN RADIANS 
Y LABEL : PROBABILITY G8 KELL 
POPULATION NO. OF se eeeee PERCENTILES=------- 
NUMBER POINTS Y MEAN ~Y Sib pEV 10-2> Ors, ee: Oe 7 
0 Log 0.16605 O-3i466 0 2.10325 04 O02 2si6 ie 
O20005 100 O.17568 0222455 @) 5.7949E _6 0.27246 
OU 2007 100 0.18497 0.31487 0 6. Isles On Zee 
D200 5s 100 O70 O2Z2776 @) 1.9134E_5 0.23033 
O.00Z 100 0.16013 0.28406 0 Si jicysieieid (6 0.13654 
POOLED STANDARD DEVIATION ESTIMATE: 0.29559 
3. RANGE EVALUATION 
TABLE DS 
RANGE EVALUATION 
PARAMETER VALUES: (100);(700);(1300);(1900) 
POPULATIONS >: RG14;RG13;RG12;RG11 
WEIGHTS : O 
% EREEL >: RANGE IN METERS 
Y ECE : PROBABILITY OF KILL 
POPULATION NO. OF = ~~ 6speeeeccece PERCENTILES**22e2e2+-%-%--- 
NUMBER POINTS Y MEAN Y SIDS EV 0.225 aS) Os7s 
100 16 O.e06l3 O31 907) 0.72854 8.6001E-1 9.4777E-1 
700 20 O29 0223669 0.021002 1.3836E=-1 3.22332 —1 
1300 50 0.046666 0.15413 0 2.0463E-9 2.12226-2 
1900 50 0.0094727 0.065348 0 0.Q000E0 6. 1062-16 


POOLED STANDARD DEVIATION ESTIMATE: 0.16975 


By 


4. EFFECTS OF TARGET MANEUVERS 


a. Maneuvering Target 


PARAMETER VALUES : 


POPULATIONS 
WEIGHTS 
A” LABEL 
Y LABEL 
FOrULATION NO. OF 
NUMBER POUNTS 
100 50 
700 20 
1300 50 
1900 50 


TLABIEE 6 
MANEUVERING TARGET 


(100) ;(700); fry L300) 
RG14;RG13;RG12;RG1 


: RANGE IN METERS 
; PROBABILITY OF KILL 


PeebED STANDARD DEVIATION ESTIMATE: 0.17357 


b. Non-maneuvering Target 


PARAMETER VALUES: 
POPULATIONS 
HELGHTS 

xX LABEL 

PecABEL 


POPULATION NO. OF 
NUMBER 


100 50 
700 50 
1300 50 
1900 a0 


POrNis «1 MEAN 


: O 


=<<S-22->— PERGHNT LILES ==—so—a----"- 
Y MEAN HC SB) Ba.) ON 8, Oo O75 
0.80613 eo 0.72854 S. O00TE =) SoS ITA 
O23 Or2 5569 O.021062 Ee ey lels— | Beco e- 1 
O05) 466 0.15441 0 3) AE ap fo oUSk=3 
0.0094727 0.065348 0 0.Q0000E0 Se OoZk- 16 
TABLE 7 
NON-MANEUVERING TARGET 
(TGC (700) sel 300 771900) 
NOM1 ; NOM2 ; NOM3 ; >NOM4 
: RANGE IN METRES 
SP ROBAB VEL] WeOr shine 
=“<s- S525 PERCENLILES= ---<s22- 
eo LD) DEV Oe: Ons On 
0.97438 0.009994 Ze G7 O2E=. 0. 97008 O.go113 
0.46539 0225852 Me OOSE =) 0.41486 0.6943 
Oro. 26 0.27431 § 307 1E-11 G. O013 258.9 0. 2ene1 
0.012354 0.075884 0.0000E0 0 


BeOOlED STANDARD DEVIATION ESTIMATE: 0.19231 


=p 


>; EFFECTS OF RADAR ERROR 


TABEESS 
RADAR ERRNO 


PARAMETER oe (0) -( 50) (0G) aie or 
POPUEATIONS RD21;RD22;RD23;RD24 
WEIGHTS 1, 
Label : PERCENT RADAR ERROR 
Y LABEL : PROBABILITY Vora Er 
POPULATION NO. OF 4. °.}»§ (320) 66-6 PERCENTILES=--=="=—2m 
NUMBER POINTS - Y MEAN Y "S22 SDEV 02.25 OF O27s 
0 60 0.3078 0.41684 1.497 9E=12 ool aye 0.804 
50 60 0.2246 0.33007 0.Q0000E0 0.0000149 0.44578 
100 60 0. T3527 0.25246 1.3878E-17 0.00014144 0.12566 
150 60 0.083055°70 2725 0. OOU0EG 0.0017942 0.07882m 
POOLED STANDARD DEVIATION ESTIMATE: 0.29038 
6. EFFECTS OF MUZZLE VELOCITY AND RATE OF FIRE 
a. Nfuzzle Velocity 
TAB IEE 9 
MUZZtE VELOC 
PARAMETER rane (1000) ;(1500); AAO) CZ 500) 
FOPULATIONS MVi- MV2; Mya; 
WE DGHES ee, 
nh LABEL : MUZZLE VELOCITY (METERS/SEC) 
¥ “LABEL : PROBABILITY OF VKIEE 
POPULATION NO.OF we.) | ee eee - PERCENT ILES <= --—— ue 
NUMBER POINTS .¥ MEAN ¥ SID GeV 0325 Oe) 0.75 
1000 60 O,11223) 0.23477 0 0.000029386 0.055264 
1300 60 0.24193 20. 29665 0.0022309 .UmiZ2 0.393e2 
2000 60 0.35197 0.27194 O2.09949 O751493 0.57484 
2500 60 0.51254 "O22572 On28252 0.5894 0. 703i 


POOLED STANDARD DEVIATION ESTIMATE: 0.26596 





b. Rate of Fire 


Waele 200 
Ree ler r rine 


PARAMETER VALVES: (400) ;(600);(800) ; pee 


POPULATIONS Paz FR22; FR23;FR24 

WEIGHTS a8) 

A LABEL : FIRING RATE(RDS/MIN) 

Pe eeSEL PRODAS id tyr eK LL 

POPULATION NO. OF ae al PERGENDLERS == ca=<-=—=< 

NUMBER POC Ss awe an 51D DEY Ox2> OR Ons 
400 60 Qelols Deo o13 O[OCtseCrs O07 55 et9513 
600 60 OG. 21 eZee 2/047 0.00079402 0.058486 0.39343 
800 60 O,30b9cme 50679 0.0028055 Weascol Os See eagle 
1000 60 0.40055 0.32489 ©7033221 0.38058 G274919 


meOlED STANDARD DEVIATION ESTIMATE: 0.29132 


ty 
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